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Statistical properties of cross sections are studied for an open system of interacting fermions. The 
description is based on the effective non-Hermitian Hamihonian that accounts for the existence of 
open decay channels preserving the unitarity of the scattering matrix. The intrinsic interaction is 
modelled by the two-body random ensemble of variable strength. In particular, the crossover region 
from isolated to overlapping resonances accompanied by the effect of the width redistribution creat- 
ing super-radiant and trapped states is studied in detail. The important observables, such as average 
cross section, its fluctuations, autocorrelation functions of the cross section and scattering matrix, 
are very sensitive to the coupling of the intrinsic states to the continuum around the crossover. A 
detailed comparison is made of our results with standard predictions of statistical theory of cross 
sections, such as the Hauser-Feshbach formula for the average cross section and Ericson theory of 
fluctuations and correlations of cross sections. Strong deviations are found in the crossover region, 
along with the dependence on intrinsic interactions and degree of chaos inside the system. 
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I. INTRODUCTION 

Information on properties of quantum mesoscopic sys- 
tems comes mostly from various reactions where the sys- 
tem plays a role of a target. At high density of intrinsic 
states, the dynamics of realistic systems of interacting 
constituents becomes chaotic. Onset of chaos immensely 
complicates the details of the scattering process as re- 
flected in the paradigm of compound nucleus. At low 
energies, the long-lived resonance states are exceedingly 
complex superpositions in the basis of independent parti- 
cles. As energy increases, the scattering pattern evolves 
from the set of narrow isolated resonances to overlapping 
resonances and strongly fluctuating cross sections. Since 
the individual properties of resonances cannot be pre- 
dicted, only statistical description is practical and sen- 
sible. The average cross sections are usually described 
according to Hauser-Feshbach [l|, while the fluctuations 
and correlations of cross sections are treated in terms of 
Ericson theory [i,|l,l3- 

Standard theory of statistical reactions does not an- 
swer the question of interplay between reactions and in- 
ternal structure determined by the character of interac- 
tions between the constituents. Here more detailed con- 
siderations are required based on the generalization of 
the shell model of nuclear reactions Such an exten- 
sion introduces statistical assumptions concerningintrin- 
sic dynamics and its coupling to the continuum [6|,|7| with 
the extreme limiting case of intrinsic chaos modelled by 
the Gaussian Orthogonal Ensemble (GOE) [sl]. An im- 
portant element is still missing here, namely the transi- 
tional regime from separated to overlapping resonances. 
The consistent description based on the continuum shell 



model ra, as well as more phenomenological approaches 
H, El , indicate the presence of a sharp restructur- 
ing of the system when the widths of resonances become 
comparable to their energy spacings. This phenomenon 
carries features of a quantum phase transition with the 
strength of continuum coupling playing the role of a con- 
trol parameter. 

As was clearly observed in the shell model framework 
[T^ l , the distribution of resonance widths rapidly changes 
in the transitional region in such a way that a number 
of very broad resonances equal to a number of open de- 
cay channels absorb the lion's share of the total width 
of all overlapped resonances, while the remaining states 
become very narrow. The corresponding theory was sug- 
gested in [l^ [3, [il], where the mechanism of this re- 
structuring was understood to be associated with the na- 
ture of the effective non-Hermitian Hamiltonian [5;] that 
describes the intrinsic dynamics after eliminating the 
channel variables. The factorized structure of this Hamil- 
tonian, in turn, is dictated by the unitarity of the scat- 
tering matrix [l6| . One can coinpare this phenomenon to 
the classical factorized model [131 of a giant resonance, 
where the collective strength of many particle-hole states 
is shifted in energy and concentrated at a speciflc com- 
bination of excited states. In spite of formal analogy, 
physics under study here is different. The concentra- 
tion of widths on a few broad states can be described as 
collectivization along the imaginary axis in the complex 
energy plane. The driving force of this restructuring is 
the presence of open decay channels and interaction of in- 
trinsic states through the continuum. The intrinsic Her- 
mitian interaction is present as well and should be fully 
accounted for; one of the goals of our study is to under- 
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stand the dependence of the continuum picture on the 
strength and character of interactions inside the closed 
system. The interplay of two collectivities is an interest- 
ing subject [H, [31 practically important in relation to 
the so-called pygmy-resonances in loosely bound systems 

The segregation of short-lived broad resonances from 
long-lived trapped states was shown to be similar to the 
superradiance [2l| in quantum optics induced by the cou- 
pling of atomic radiators through the common radiation 
field, an analog of coherent coupling of many overlapped 
intrinsic states through continuum. Later a general char- 
acter of the phenomenon was demonstrated for systems 
with GOE intrinsic dynamics and many open channels 
[23 , [SS*! . Modern versions of the shell model in contin- 
uum [24,, 25] are based on the effective Hamiltonian and 
naturally reveal the superradiance phenomenon as an im- 
portant element. The transition to this regime should be 
taken into account in all cases when a physical system is 
strongly coupled to the continuum, see for example (26l | 
and references therein. 

The segregation of scales is srcctacularly seen in level 
and width statistics [13, [H, [g^, HI]. Even for GOE 
intrinsic dynamics, the probability P{s —^ 0) of very 
small spacings between the centroids of resonances does 
not vanish because of energy uncertainty of unstable 
states. The width distribution reveals the separation of 
the "cloud" of superradiant states far from the real en- 
ergy axis, while the trapped states are clearly accumu- 
lated near the real axis. Some features of the picture 
are however sensitive to the character and strength of 
intrinsic interactions as we have demonstrated recently 
[2^ . The analysis included the two-body random ensem- 
ble of variable interaction strength in a Fermi system of 
shell-model type and the GOE as an extreme limit of 
many-body random interaction. 

In the present work, using the same framework as in 
[29! , we study the interplay between the intrinsic dynam- 
ics and statistical properties of cross sections comparing 
the results with those of conventional approaches, namely 
Hauser-Feshbach average cross sections and Ericson fluc- 
tuations and correlations. In particular, we show that 
the assumption that fluctuations of the resonance widths 
are negligible for a large number of channels is not cor- 
rect. We also show that the elastic enhancement fac- 
tor strongly depends on the degree of chaoticity inside 
the system, thus leading to deviations from the Hauser- 
Feshbach formula. 



II. THE MODEL 

A. Hamiltonian 

We consider a system of n interacting fermions on m 
mean- field orbitals. A large number, N = TO!/[n!(m — 
n)!], of intrinsic many-body states \i) comprise our 
Hilbert space. In our simulations we take n = 6, to = 12 



that provides a sufficiently large dimension N = 924. 
The states are unstable being coupled to M open decay 
channels. The dynamics of the whole system is governed 
by an effective non-Hermitian Hamiltonian [ij, [2^ 
given by a sum of two N x N matrices, 

M 

n = H-'-W; W^.^Y.'^^- (1) 

Here and below the intrinsic many-body states are la- 
belled as i^j,.-. and decay channels as a,b,c.... In Eq. 
(H)), H describes Hermitian internal dynamics that in re- 
ality also can be influenced by the presence of the con- 
tinuum [23, [2^, [s^l , while is a sum of terms factorized 
in amplitudes A1 coupling intrinsic states \i) to the chan- 
nels c. Under time-reversal invariance, these amplitudes 
can be taken as real quantities so that both H and W 
are real symmetric matrices. 

We model H by the two-body random ensemble 
(TBRE) assuming the intrinsic Hamiltonian H in the 
form H = Hq + V, where Hq describes the mean-field 
single-particle levels \v), and ^ is a random two-body in- 
teraction between the particles [3l|. The single-particle 
energies, e,y, are assumed to have a Poissonian distribu- 
tion of spacings, with the mean level density l/do- The 
interaction V is characterized by the variance of the two- 
body random matrix elements, (V^^ ,^^.^^ ,^^) = Vq. With 
no interaction, vq — 0, the many-body states have also 
the Poissonian spacing distribution P{s). In the oppo- 
site extreme limit, do = 0, corresponding to infinitely 
strong interaction, A = vo/do — > 00, the function P{s) 
is close to the Wigner-Dyson (WD) distribution typical 
for a chaotic system Q. Following Ref. the critical 
interaction for the onset of strong chaos can be estimated 
as 

_ fcr _ 2(to - n) 

where Ng — n{m — n) -I- n{n — 1)(to — n)[m — n — l)/4 is 
the number of directly coupled many-body states in any 
row of the matrix Hij. Thus, we have Acr ~ 1/20, and 
often we perform the simulations with the value A = 1 /30 
slightly lower than Acr- In parallel, we also consider the 
intrinsic Hamiltonian H belonging to the GOE that cor- 
responds to a many-body interaction, when the matrix 
elements are Gaussian random variables, (Hfj) ~ 1/N 
for i ^ j and (Hf^) = 2/7V for i = j. 

The real amplitudes A'^ are assumed to be random in- 
dependent Gaussian variables with zero mean and vari- 
ance 

{AIA'^') = S,,6^^' ^. (3) 

The parameters 7*^ with dimension of energy character- 
ize the total coupling of all states to the channel c. The 
normalization used in Eq. ([3]) is convenient if the energy 
interval ND covered by decaying states is finite. Here D 
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is the distance between the many-body states in the mid- 
dle of the spectrum, D = l/p(0), where p{E) is the level 
density, and E = Q corresponds to the center of the spec- 
trum. We neglect a possible explicit energy dependence 
of the amplitudes that is important near thresholds and 
is taken into account in realistic shell model calculations 
[25I [30| . The ratio ^'^ /ND characterizes the degree of 
overlap of the resonances in the channel c. We define the 
corresponding control parameter as 



7r7 

2ND' 



(4) 



The transitional region corresponds to k'^ ~ 1. Varying 
the intrinsic interaction and, therefore, the level density 
p, we renormalize correspondingly the absolute magni- 
tude of the widths, 7, in order to keep the coupling to 
continuum given by Eq. ([4]) fixed. 



B. Scattering matrix 

The effective Hamiltonian allows one to study the cross 
sections for possible reactions 6 — > a. 



a'^iE) = \T'"'{E)Y 



(5) 



(our cross sections are dimensionless since we omit the 
common factor 7r/fc^). In what follows we study both 
the elastic, b — a, and non-elastic, b ^ a, cross sections. 
Ignoring the smooth potential phases irrelevant for our 
purposes we express the scattering amplitude of the re- 
action, T*"*, in terms of the amplitudes A?, 



T'"'{E) 



1 



N 

^ ' \E-n, „ 



4. 



(6) 



Here the denominator contains the total effective Hamil- 
tonian ID) including in this way the continuum coupling 
W to all orders. 

We can also write T^"-{E) in a different way, diago- 
nalizing the effective non-Hermitian Hamiltonian 7i. Its 
eigenf unctions |r) and (f| form a bi-orthogonal complete 
set. 



n\r)^£r\r), (f|H=(f|f;, 
and its eigenvalues are complex energies. 



(7) 



(8) 



corresponding to the resonances with centroids Er and 
widths r^. The decay amplitudes A\ are transformed 
according to 



(9) 



and the transition amplitudes are given by 



T'^'^E) = 



N 



(10) 



The bi-orthogonality of the transformation ensures that 
the statistical properties ([3]) of the ensemble of the am- 
plitudes are preserved. 



~ h hi 



(11) 



Introducing the matrix in channel space analogous to 
what is routinely used in the resonance data analysis. 



1 1 7. 
^ ^ 2^E-H^ ' 



(12) 



where the denominator includes only intrinsic dynamics 
and A is the N x M matrix of the transition amplitudes 
A'j, one can relate it to the transition amplitude T and 
the scattering matrix. 



in the explicitly unitary form 



1-iK 
1 + iK' 



(13) 



(14) 



The poles of the if-matrix are real eigenvalues Ea of the 
intrinsic Hcrmitian Hamiltonian H. The complex eigen- 
values ([8]) coincide with the poles of the S'-matrix and, 
for small values of 7, determine energies and widths of 
separated resonances. With an increase of 7, the reso- 
nances start to overlap leading to specific features of the 
scattering process which are of our main interest. 



III. AVERAGE SCATTERING MATRIX 



The scattering matrix (|T4|) averaged over the ensemble 
of decay amplitudes ^ with accuracy of 1/A^, is given 
by iU] 



{S) = 



1 - i{K) 
l + i{K)' 



(15) 



Following Refs. [iJjlS^], we assume, in concordance with 
the statistical ansatz that the decay amplitudes and 
the eigenvalues E^ of the intrinsic Hamiltonian are sta- 
tistically independent. This assumption is satisfied for 
the GOE and Poissonian ensembles, and here we assume 
that this is also true for the TBRE. 

Since the mean field many-body basis \i) and the 
eigenbasis |a) of the intrinsic Hermitian Hamiltonian are 
related by a real orthogonal transformation, the anti- 
Hermitian part W of the effective Hamiltonian in the 
basis I a) still has a factorized form in terms of the new 
amplitudes B^; similarly to Eq. (dU, the correlation 
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function of the amplitudes B coincides with that given 
by Eq. ([3]). Thus, the average matrix determined by 
the statistical properties of the decay amplitudes reduces 
to 



{K-\E)) = -y 



^ ^ab T_ 



1 



2N ^ E-Ea 



(16) 

The last sum in Eq. (fTB|) is the trace of the intrin- 
sic Green function \/{E — H). For energy E inside 
the spectrum of H we should understand it as a lim- 
iting value, E ^ E + iO. The trace of the imaginary 
part of the Green function determines the level density 
p{E) = J2a ^(-^ ~ ^oi) for the Hamiltonian H, and 



E 



1 



E-E^ 



iO 



1 



E- Er 



-mp{E). (17) 



The principal value part is a smooth function of energy 
that vanishes in the middle of the spectrum; as a result, 
in this vicinity 

{K^'^) = -inS'^" ^ piO) = -iS'^^K'^, (18) 

where k"' is defined by Eq. fl]). 

The average scattering matrix p3|) takes the form 



1 - k'^ 

H- K° 



(19) 



that depends only on our parameter of continuum cou- 
pling independently of the intrinsic interaction strength 
A. This independence is due to the fact that we normal- 
ized K to the mean level spacing, D, but it is important 
to stress that in the shell model D does depend on the 
intrinsic interaction [3^. Eq. fTO]) implies that the trans- 
mission coefficient. 



T" = 1- 
is also independent of this strength, 



(1 -f K'^y 



(20) 



(21) 



The transmission coefficient in the channel a, T", is max- 
imum (equal to 1) at the critical point of this channel, 
k" = 1, when the average S'-matrix vanishes. Thus, 
K = 1 determines the so-called perfect coupling regime. 
We will study the statistical properties of cross sections 
as a function of the intrinsic interaction strength A and 
continuum coupling parameter k, both below the critical 
point, K < 1, and after the superradiance transition has 
occurred, k > 1. 



TBRE is Gaussian for large enough particle number, n, 
and orbital number, m, [8|; its width depends on A. In 
order to compare different ensembles, we restrict our sta- 
tistical analysis to a small energy interval with a constant 
level density at the center of the real spectrum of the 
complex eigenvalues of Ti.. This interval should be small 
enough in order to neglect the energy-dependent differ- 
ence of the density of states among the ensembles, but 
large enough with respect to the widths in order to con- 
tain a statistically meaningful number of resonances. 

For a model with a finite resonance number, it is im- 
portant to avoid edge efi'ects, see discussion in [33 |. 
The energy interval subject to statistical analysis should 
be also at a distance of at least several widths away from 
the edges. A rough estimate l^o'] goes as follows: for M 
equivalent channels, T/D oc M, and the distance from 
the center to the edges is ND/2, then ND/2 > T/D that 
implies M/N ^ 1/2. This shows that the ratio of the 
number of channels to that of resonances must be small 
in order for the results to be model independent. With 
this choice, the model will be essentially equivalent to an 
infinite resonance model with a constant level density, 
apart from a narrow interval around the critical value. 
For K = 1, with an infinite resonance number, the aver- 
age widths should logarithmically diverge, in agreement 
with the Moldauer-Simonius expression 11]. In our finite 
model, the results become model-dependent in a narrow 
interval near k = 1. Our approach is still appropriate 
for a comparison with the predictions of Ericson fluctua- 
tion theory derived for an infinite resonance model with 
a constant level density. 

The results of numerical simulations presented below 
refer to the case of TV = 924 internal states and M 
equiprobable channels, = k. The maximum value of 
M we considered is M = 25 so that M/N fa 2- 10"^. For 
any value of k, we have used a large number of realiza- 
tions of the Hamiltonian matrices, with further averaging 
over energy. 



V. ERICSON FLUCTUATIONS 

The starting point of the conventional theory 0, H, 01 
can be summarized as follows. The scattering amplitude 
T''''{E) = (T"''(£;)) + T^'^iE) is divided into two parts, 
an average one, {T°-^{E)), and a fluctuating one, T^^{E), 
with 



{Tt{E)) ^ 0. 



(22) 



Note that in our statistical model we have (Tinci) for 
inelastic channels, while 



IV. COMPARING THE ENSEMBLES 

The density of states of the GOE ensemble follows the 
famous semicircle law, while the density of states of the 



(Tel) = -z(l - (5)) = -2i 



(23) 



for clastic channels. 

With statistical independence of poles (resonance en- 
ergies) and residues (resonance amplitudes), Zr = A'^A^, 
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we obtain = (z^) + 5zr^ so that 



{Zr) + 5z 



E - Er + iTr/2' 



(24) 



In the regime of overlapping resonances, (F) > D, and 
assuming ah widths of the same order, Tr ~ (F), the 
average part can be computed similarly to Eq. (I17p , sub- 
stituting the sum by the integral, 



p{Er){Zr)dEr .T^{Zr) 

E - Er + i{T) /2 ~* D 



(25) 



where a constant level density, p{Er) = 1/-D, is assumed. 

The average cross section, cr = |Tp, also can be di- 
vided into two contributions, 

(a) = (|Tp) = |(T)p + (|Tflp). (26) 

The two terms in Eq. (pS)) are interpreted as 



(O'dir, 



(27) 



where the direct reaction cross section, (o'dir), is deter- 
mined by the average scattering amplitude only, while 
(cTfl) is the fluctuational cross section (also called in the 
literature the compound nucleus cross section) that is 
determined by the fluctuational scattering matrix. 

For overlapping resonances, (F) > D, the following 
conclusions were derived concerning the scattering am- 
plitude and the statistical properties of the cross sections. 

(A). The average fluctuational cross section As- 
suming that Fr « (F) for a large number of channels, i.e. 
fluctuations of the widths around their average value are 
small, Var(F)/(F)^ ^ 1, the average fluctuational cross 
section, (cta) ~ (TflTg), can be written as 



E 



Sz*SZr' 



-{E-Er+ t/2{T)){E -E^- i/2{T)) 



(28) 

where the substitution F^ ~ (F) was used. Now the 
averaging over energy is applied. 



{F{E))^^ I dEF{E). 



The integration leads to 
, , 27ri 



AE ^ Er' - Er + i(T)' 



(29) 



(30) 



Now we assume that dz^ are uncorrelated random quan- 
tities with the statistics independent of r, {Sz*dzr') = 
6rr'{\Sz\'^). The absence of correlations between the am- 
plitudes Szr gives 

2n{\6z\^) AE 



(B) . Variance of the cross section, Va.r{a) = (cr^) — (cr)^. 
The derivation can be performed under more general as- 
sumptions [3^ than those used for the analysis of average 
cross section. If the quantities 5zr/{E — £r) are indepen- 
dent complex variables, then T is Gaussian distributed, 
that is T = ^ + irj, where both ^ and 77 arc Gaussian 
random variables with zero mean. This is due to the fact 
that for (F) ^ D both ^ and 77 are the sums of a large 
number of random variables. In the conventional theory 
it is also assumed that ^ and rj have equal variance. 

Then for the fluctuating cross section we have 

(ai) = (|Tfl|4) = {T^TiT^Ti) ^ 2{a^)\ (32) 
In a more general case when (T) 7^ 0, 

Var(a) = (afl)(2(adir) + (fTfl)). (33) 

(C) . The correlation function of the scattering ampli- 
tudes is defined as 

c{e) = {T{E + e)T*{E))-\{T{E))\'' = {%{E + e)T^ {E)) . 

(34) 

Evaluating c(e) under the same assumptions as for the 
average cross sections, one obtains. 



c(e) = (cTfl) 



(35) 



e + z(F)- 

(D). The cross section correlation function is defined 



as 



C{e) = {a{E)a{E + e))-{a{E))' 



(36) 



Taking into account the Gaussian form of distribution for 
T and Eq. ipS)) . one obtains that: 

(a) the normalized autocorrelation function of cross 
sections satisfies the relation, 



Cie) |c(6)|^ 



C(0) |c(0)P' 
(b) the correlation function has a Lorentzian form, 

Cje) _ P 
C{0) P + e2' 



(37) 



(38) 



where the correlation length, I, is equal to the average 
width, 



Z=(F). 



(39) 



In the following we compare the predictions (B)-(D) 
of the conventional theory of Ericson fluctuations with 
our numerical results, paying special attention to the de- 
pendence on the intrinsic interaction strength, A. As for 
the part (A), since there are no predictions for the quan- 
tity |(5zp determining the average fluctuational cross sec- 
tion pip , the comparison will be done with the Hauser- 
Feshbach, theory widely used in the literature. 
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VI. AVERAGE CROSS SECTION 

In this section we study how total and partial cross 
sections depend on the continuum coupling, k, and in- 
trinsic interaction, A. For any value of k, we have used 
Nr — 30 realizations of the Hamiltonian matrices. For 
each realization we took into account only the interval 
[-0.2,0.2] of real energy at the center of the spectrum. 

It follows from Eq. that the average total cross 

section defined by the optical theorem. 



(atot) =2(1-Rc(5)) = - 



4k 



(40) 



depends only on the average scattering matrix and there- 
fore is independent of A and M . Since a^x = txtot for 
M = 1, the average elastic cross section is also inde- 
pendent of A for the case of one channel. The situation 
changes as we increase the number of channels. 

In order to analyze the average elastic and inelastic 
cross section, we single out the average scattering matrix 
elements in the standard form, 5°^ = (5°^) + where 
{S"^) = S^^iS"") and {Sf') = 0. The average inelastic 
cross section, a 7^ 6, is given by 



^ab|2\ 



ab\2\ 



The average elastic cross section can be written as 

(cr"'^) |1 - (S"^'^>P + 



(41) 



(42) 



Following the literature we will call (|S'g''p) the fluctu- 
ational cross section, {cr'^). For M equivalent channels 
the fluctuational cross section can be expressed with the 
use of the elastic enhancement factor. 



F 



-aa\ 



(43) 



where b ^ a. Indeed, in the case of equal channels, using 
the relation. 



we obtain 



{M - ly"' + \1 - S' 



^aa |2 



1 - |(S"''')|2 _ T 

F + M -1 ^ F + AI-1' 



(44) 



(45) 



a detailed analysis of the elastic enhancement factor for a 
separate study, here we point out that F also depends on 
K. Specifically, with an increase of k from zero, the value 
of F decreases, being confined by the interval between 3 
and 2. 

For the fluctuational inelastic cross section, with an 
increase of the number of channels the dependence on 
the interaction strength A disappears, see Fig. [1] (upper 
panel). This is in agreement with the fact that in the 
limit of large M we have {a^'') T/M, independently 
of A, see Eq. ([15)) . In contrast, the fluctuational elastic 
cross section manifests a clear dependence on A, see Fig.[T] 
(lower panel). Thus, one can directly relate the value of 
the enhancement factor F to the strength A of interaction 
between the particles. As one can see, the more regular 
is the intrinsic motion, the higher is the average cross 
section. Note that for a large number of channels the 
A— dependence of the elastic cross section is in agreement 
with the estimate, (o-g") FT/M, see Eq. (gH). 




FIG. 1: (Color online) Fluctuational inelastic (upper panel) 
and fluctuational elastic (lower panel) cross sections as a func- 
tion of K. Circles refer to the GOE case, pluses to A = 0, 
crosses to A = 1/30, and squares to A 00. The solid 
curves (top-down) correspond to the HP formula (|49|) with 
F = 2.5, 2.25, 2.0 for A = 0, 1/30 and the GOE, respectively. 
The values 2.5 and 2.25 are found numerically using the def- 
inition (l43ll. 



where T is the transmission coefficient defined in Eq. 
(ED), and 



ab\ 



FT 



F + M -r 



a^b. 



(46) 



Since the transmission coefficient T does not depend 
on A, the only dependence on A in Eqs. (|45)) and ([46]) 
is contained in the elastic enhancement factor F. The 
same seems to be correct even when the channels are non- 
equivalent, according to the results of Ref. ^3^. Leaving 



Now we compare our numerical results with the 
Hauser-Feshbach (HF) formula derived under quite gen- 
eral assumptions for both a = b and a ^ b (see Ref. Q 
and references therein). 



«) = (l + 5,,) 



(47) 



Here the last expression corresponds to our case of equiv- 
alent channels, T'^ = T. As one can see, the HF formula 
predicts F = 2, independently of the interaction strength 
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between the particles. This formula was also derived in 
in the overlapping regime for the TBRE ensemble with 
the infinite interaction, and in [s^l for the GOE ensemble. 
For finite number of channels in the limit \/M ^ 1, the 
corrected HF formula was derived in Q. For equivalent 
channels it reads, 



«) = (l + <5,,) 



T 



1 



1 

M 



(48) 



Our data confirm that for the fluctuational inelastic 
cross section the HF formula gives correct results for all 
values of A in the case of large number of channels. The 
specific case of small number of channels, for which the 
HF is not valid, will be discussed elsewhere. 

On the other hand, for the fluctuational elastic cross 
section, our data show that the HF formula works only in 
the GOE case and in the limit A ^ oo, see Fig. [1] (lower 
panel) . At finite values of A clear deviations are seen. In 
order to describe the data, we modified the HF formula 
taking into account that the elastic enhancement factor 
varies with A, 



[l + nF-l)]J(l-^) 



(49) 



As one can see, this expression gives a satisfactory de- 
scription of the data, with the numerically computed 
values of F. The problem of an analytical dependence 
of F on the interaction strength A remains open. To 
shed light on this problem, we performed a specific study 
of the elastic cross section in dependence on A for fixed 
value K = 0.8 in the overlapping regime, see Fig. 2. As 
one see, there is a sharp decrease of the cross section in 
the transition from regular to chaotic intrinsic motion, 
A ~ Aci-. This result is quite instructive since it shows 
how the scattering properties are influenced by the onset 
of chaos in an internal dynamics. The non-trivial point 
is that the analytical estimate of Acr was obtained for a 
closed system, k — 0. However, even in the regime of a 
strong coupling to the continuum, k = 0.8, this estimate 
gives a correct value for the interaction strength at which 
a drastic change of scattering properties occurs. 

VII. FLUCTUATIONS OF WIDTHS AND 
RESONANCE AMPLITUDES 

Here we discuss the conventional assumption that for 
a large number of channels the deviations of the widths 
from their average are small, Var(r)/(r)^ <C 1, and, 
therefore, for analytical estimates one can set F^ ~ (F), 
see Eq. (PSI) . It is usually said in justification of this as- 
sumption [2] that in the overlapping regime the width can 

be presented as a sum of partial widths, F^ = X]c=i ^r- 
Assuming that individual partial widths obey the Potter- 
Thomas distribution, the total width is expected to have 
a Xm distribution, so that Var(r)/(F)2 = 2/M is smafl 
for M ^ 1. In fact, it is sufficient to accept that the 
partial widths are independent random variables; then 
Var(F) cx M and (F) cx M, so that Var(F)/(F)2 cx Af-^. 



.\o.23 
D 

V 0.22 



0.21 
0.2 
0.19 
0.18 
0.17 



1,1,1,1 

I 1 
• U A 

T 1 cr 


1 1 1 , 1 1 1 1 


1^0.8 


-GOE i 

\ 1 ^ 











0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

X 

FIG. 2: (Color online) Fluctuational elastic cross section as a 
function of the interaction strength A for M — 10 and k — 0.8 
(connected circles). The horizontal line refers to the GOE 
value, and dashed vertical line shows the critical value Acr for 
the transition to chaos in the TBRE, see Eq. ([2]). 



However, in our previous work |29l | we have showed 
that for large values of k the distribution of the widths 
strongly differs from the Xm distribution. Our new 
data in Fig. 3 give more details concerning this problem. 
These data were obtained for a large number = 100 
realizations of the Hamiltonian matrices, in order to have 
reliable results. 

The data show that as k increases the normalized vari- 
ance, Var(F)/(F)^, also increases, remaining very large 
even for M — 20. Moreover, the deviations from the ex- 
pected 1 /M behavior are clearly seen signaling the pres- 
ence of correlations in the partial widths. From Fig. 3 
one can also understand how the value of Var(F)/(F)^ 
depends on the degree of intrinsic chaos determined by 
the parameter A. Specifically, for small k there is no de- 
pendence on A and Var(F)/ (F)^ decreases as 2/M for all 
the ensembles, as expected. However, as k grows, the 
dependence on A emerges: the weaker the intrinsic chaos 
(and, consequently, the more ordered is the intrinsic spec- 
trum) the larger are the width fluctuations. 

In Ref. m the widths distribution for the GOE en- 
semble was found analytically for any number of channels 
in the limit of oo and M fixed. The general ana- 

lytical result is given in terms of a complicated threefold 
integral. A simpler expression is obtained for a specific 
case M = 2, for which the distribution of the widths is 
expressed as a double integral. 



dv 



1 



(50) 

(2/r) - 1 and 



where C = l/{2./f~^), y = ttF/D, g 
T is the transmission coefficient, so that from Eq. ((2T|) 
we have 5 = (1-1- k^)/2k, or k 



g ± _ I Note that 
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FIG. 3: (Color online) Normalized variance of the width as a 
function of the number of channels M, for different coupling 
strengths k (connected symbols are the same as in Fig. [TJ. 
While for small coupling, k = 0.01, the variance decreases 
with the number of channels very fast in accordance with the 
expected x^-distribution (dashed line), for large couplings. 



K = 0.5 and 0.9, 
dependence. 



the behavior is different from the 1/M- 



derived from Eq. ([50]) for this region has the form, 

2 



Var(r) 



2(2 + 7r) 



v/2(5 + l)(5-l) V 9 + 



In 



.9-1 



1 



1, 



(51) 



and is shown in Fig. 4 by the dashed curve. From the 
above relation for M = 2 one can obtain that the nor- 
malized variance diverges as 



Var(r) 



oc 



(1 



^[\n{l^^)]\ (52) 



Our numerical simulations confirm that the divergence 
remains for any number of channels and for any value of 
A. Thus, contrary to the traditional belief, the variance 
of widths does not become small for a large number of 
channels. Another result is that the assumption of the 
absence of correlations between the resonance amplitudes 
Szr and the widths F^ seems to be incorrect in the region 
of a strong resonance overlap. Indeed, the data reported 
in Fig. 5 demonstrate that in contrast with the case of 
weak coupling, k = 0.001, for a strong coupling there 
are systematic correlations between Szr and F^. These 
correlations are increasing with an increase of the cou- 



v + g — 2fi > and for g = 1 we have k = 1. The result 
of numerical integration in Eq. (I50p is shown in Fig. [4] 
by a solid curve. The comparison with the normalized 
variance of the widths for the GOE case and M = 2 
(circles) shows a good agreement except for the vicinity 
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FIG. 5: (Color online) Absolute squares of resonance ampli- 
tudes, |5zr-|'^, versus the widths Tr for the GOE with M — 20. 
As K increases, the correlations between \Szr\'^ and grow. 



FIG. 4: (Color online) Numerical data for the normalized vari- 
ance of the widths vs k for GOE and M = 2 (circles) , in com- 
parison with the result of numerical integration of Eq. (|50p 
(solid curve), and with Eq. (|5ip (dashed curve), see in the 
text. 

A specific interest is in the behavior of the variance at 
the transition region, k « 1. The analytical expression. 



VIII. STATISTICS OF CROSS SECTIONS 

A. Distribution of fluctuational cross sections 

According to the standard Ericson theory, the fluctuat- 
ing scattering amplitude can be written as = 7] + 
where i] and ^ are Gaussian random variables with zero 
mean and equal variances. Since the fluctuating cross 
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section is given by 



(53) 



then (jfl should have a distribution with two degrees 
of freedom, that is an exponential distribution, 



P(x) = e 



<7fl 



(54) 



This should be valid both for the elastic and inelastic 
cross sections. 
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FIG. 6: (Color online) Distribution of the inelastic fluctua- 
tional cross section for the GOE and for the A = case, for 
M = 10; 25 number of channels and fixed k = 0.9. 

In Figs. [5] and [7] we show the distribution of fluctua- 
tional cross sections for the elastic and inelastic cross sec- 
tions, with two different numbers of channels, M = 10 
and M = 25. Analyzing these data, one can draw the 
following conclusions. First, for the inelastic cross sec- 
tion the data seem to follow the predicted exponential 
distribution. It should be noted, however, that a more 
detailed analysis with the help of the x^-test reveals the 
presence of strong deviations. 

The situation with the fluctuational elastic cross sec- 
tion is different due to strong deviations from the ex- 
ponential distribution occurring even for a quite large 
M — 25. The fact that large deviations from the conven- 
tional theory (for finite values of M) should be expected 
in the elastic case were recognized also in Refs. ,38, 39]. 
The comparison between M — 10 and M = 25 cases in- 
dicate that it is natural to assume that with a further 
increase of M both the distributions will converge to the 
exponential one. It is important to note that there is a 
weak dependence on the interaction strength A between 
the particles. This is confirmed by a closer inspection 
of the data of Fig. [71 Specifically, the data clearly show 
that there is a systematic difference for the two limiting 
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FIG. 7: (Color online) The same as Fig. 6, but for the elastic 
fluctuational cross section. A clear difference from the expo- 
nential distribution is seen in the tails. 



cases of zero and infinitely large values of A. Our results 
for the normalized variance (see below), indeed, confirm 
a presence of this weak dependence on A. 

The above data for the distribution of the normalized 
fluctuational cross section may be treated as a kind of 
confirmation of the Ericson fluctuation theory. However, 
it should be stressed that if we are interested in the fluc- 
tuations of non-normalized cross sections (at least, for 
elastic cross sections), one should take into account the 
dependence on A. Currently no theory allows one to ob- 
tain the corresponding analytical results, even for the sit- 
uation where the number of channels is sufficiently large. 



B. Fluctuations 

Here we compare our results for the variance of cross 
sections with the Ericson fluctuations theory 0] , and with 
more recent results for the GOE [H, [s^. According to 
the standard predictions, the variance of fluctuations of 
both elastic and inelastic cross sections. 



Var(a'^^) = ((fi"'' - (a'^''))^). 



(55) 



is directly connected to the average cross sections by 
Eq. p3)) . It is useful to express the variance of the cross 
sections in terms of the scattering matrix. In our statis- 
tical model for a ^ 6, (T'^^) = = 0. Therefore, the 
variance of the inelastic cross section reads 



Var(cr 



ah\ 



I „ab\ 1 



a6|2\2 



(56) 



For the elastic scattering one can write {T"'"') = 
{S-'^)), so that (Tdir = |1 - and afl = (15^1') 
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Therefore, for the variance of the elastic cross sections, Comparing Eq. ([58)1 with the standard predictions, Eqs. 
one obtains, (|56l57p . one can see that they are correct if: 



Var(a'^'') = 2{ar){^Z) + K">' = (i) (l^f H - 2(|5f p)^ ^ Q, 



(57) 



(ii) (^ri^^f 



afc|2\ 



0. 



(59) 




FIG. 8: (Color online) Variance of inelastic (upper panel) and 
elastic (lower panel) cross sections, see Eq. (I27|) . for M = 10 as 
a function of k, for different interaction strengths: the GOE 
(circles), A = (pluses), A ^ oo (squares), and A = 1/30 
(crosses). For the comparison the theoretical curve for MT ^ 
1, obtained for the GOE from Eqs. (|60l6ip (solid curve) is also 
added. Note that the maximum of Var(cr'^') is shifted from 
K = 1 due to the presence of direct processes. 

The discussed above conventional predictions and our 
analysis of the average cross section in Sec. VI imply that 
the variances of cross sections depend on the intrinsic in- 
teraction strength A through the average cross sections. 
Our data for a relatively large number M = 10 of chan- 
nels, see Fig. [H indeed, correspond to this expectation. 
Since for the inelastic scattering the average cross section 
does not depend on the interaction strength, it is quite 
expected that the same occurs for the variance of the in- 
elastic cross section. Our data confirm this expectation. 
On the other hand, the variance of the elastic cross sec- 
tion reveals a clear dependence on the value of A. As one 
can see, this dependence is quite strong in the region of 
strongly overlapped resonances, for k « 1. 

Let us now compare our data with the exact expres- 
sions for the variance of cross sections. To do this, it is 
convenient to express this variance in terms of the scat- 
tering matrix, 

Var(cT''^) = {\Sf\') {\Sf\Y+ 



(i-(^""*))(i^g°|25r 



-2|l-(S'"")|2(|S'°' 



6|2\ 



These properties are consistent with the Gaussian char- 
acter of the distribution for the fluctuational scattering 
matrix. 

The analytical expressions for the variance of elastic 
and inelastic cross sections were obtained in Refs. (ssl . 
I39I for the GOE case, any number of channels and any 
coupling strength with the continuum. However, simple 
expressions were derived only for MT ^ 1. Even under 
such a condition, the analytical results show deviations 
from the conventional assumptions. Specifically, it was 
found, 



(i) (l^fl 



b|4\ 



2{\Sf\''? 



b\2 



and 



(l + 7(5"'')[6-4(r" + T'') + r2; 



(ii) (l^n'^D = -8(5"°*) 



2(T°r'') 



a^3 



(5l + l)2' 



(60) 



(61) 



(58) 



where 5i = ^ ^2 = El^")' and = (^2 + l)/(5i + 

The theoretical values for the variance of the cross sec- 
tions obtained from Eqs. (|60l6ip and from the HE for- 
mula, through Eq. ([55)) . are shown in Fig. [5] by sohd 
curve. The agreement is good for the GOE case in the 
strong coupling regime, as expected. As we can see from 
Eqs. (|60l6ip . assumptions of Eq. are valid for large 
Si . In particular it was shown in 38, that the ratio 
Var((Tfl)/ (fJfl)^, being equal to one in standard theory, sig- 
nificantly differs from unity in the range 10 < MT < 20, 
where this theory is expected to be valid. 

It is now instructive to see how the normalized variance 
of the cross sections depends on the number of channels, 
see Fig. [HI According to the Ericson prediction, the ratio 
of the variance to the square of the mean of the cross 
sections has to be 1 for strongly overlapped resonances, 
in the limit of large number of channels. As one can see, 
the data for the inelastic scattering roughly confirm this 
prediction. It is not a surprise that practically there is no 
dependence on the strength A of interaction between the 
particles. On the other hand, there is a small systematic 
difference from the predicted value, that emerges for all 
values of A, as well as for the GOE case. One can expect 
that this difference disappears for a much larger number 
of channels. 

A more interesting result arises for the elastic cross 
■ section. As one can see from Fig. [SI in this case there is 
a strong difference from the limiting value of one even in 
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FIG. 9: (Color online) Normalized variance for k = 0.9 ver- 
sus the number of channels for the elastic and inelastic cross 
sections: GOE case (connected circles), A = (connected 
pluses), A — > oo (squares), and A = 1/30 (crosses). The theo- 
retical value in the Ericson theory, Var((T)/(a)^ — 1, is shown 
as a dashed line. 



the case when the number of channels is relatively large, 
M = 10 — 20. The data clearly indicate that for the 
applicability of the standard theory one needs to have a 
very large number of channels, at least larger than M — 
25. 

Another instructive observation is a weak dependence 
of the normalized variance on the interaction strength 
A. This result is in agreement with the data reported 
in Fig. [7] for the distribution of individual values of the 
cross section, where a systematic deviation can be seen 
when comparing the GOE case with the case of A = 0. 
We would like to stress that the weak A-dependence is in 
contrast with a strong dependence occurring for the non- 
normalized variance, see Fig. [S] One can treat this effect 
as manifesting that both the variance and the square of 
the cross section average depend on A practically in the 
same way. Therefore, their ratio turns out to be almost 
independent on A. As one can see, although the standard 
predictions are not correct for the non-normalized vari- 
ance, they are in a good correspondence with the data 
for the normalized variance. 



C. Correlation functions 

Here we compare our results for the correlation func- 
tion of cross sections and the scattering matrix with the 
standard predictions, see Section V. The correlation func- 
tions for the cross section and for the scattering matrix 
were computed according to Eqs. p4l36p for the elastic 
and inelastic cross sections. The correlation lengths for 
the cross section, la-, and for the scattering matrix, Is, are 
defined as the energy for which the correlation function 
is 1/2 of its initial value. Our results can be summarized 



as follows. 

(A) For large M, we found ~ Is for any interac- 
tion strength A. On the contrary, for smaller M, our 
data show that Is < la, and this difference grows for the 
weaker interaction between the particles, A. 

(B) For large M, the correlation functions are 
Lorentzian for all A, while for a small number of channels 
the correlation function is not Lorentzian, in agreement 
with the results of [4d|. Moreover, for any M, the corre- 
lation length is different from the average width, as one 
can expect due to Eq. ((39)) . apart from the region of 
small K, see Fig. [TU] For a large number of channels, the 
correlation length is determined instead by the transmis- 
sion coefficient through the Weisskopf relation, see [s^l 
and references therein, 
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FIG. 10: (Color online) Average width, (F), and correlation 
length, Icr, normalized to the mean level spacing at the cen- 
ter of the spectra, versus k for AI — 1 (upper panel) and 
Af — 10 (lower panel). Solid curves show the MS-expression 
(|6ip . Open circles refer to the GOE case, pluses to A = 0, 
squares to A ^ cxd, and crosses to A = 1/30, all for the nor- 
malized average width. Full triangles stand for the normal- 
ized correlation length at A = 1/30. The dashed line shows 
the Weisskopf relation (|60|) . 

In Fig. 1111 it is shown that, for a large number of 
channels, the elastic correlation length is in agreement 
with Eq. (|S^ for all values of the interaction strength, A. 
The same occurs for the inelastic correlation length. 

The Weisskopf relation ([5^ has been also derived in 
Ref. for small values of the ratio m = M/N, in the 
overlapping regime for the TBRE ensemble with the infi- 
nite interaction, as well as in Ref. f?\ for the GOE ensem- 
ble. In [37] the correlation function for the GOE ensem- 
ble was computed also when m = M/N is not small, and 
the deviations from Eq. (j62p and from the Lorentzian 
form of the correlation function were found. This is not 
in contrast with our results for the case of small ratio of 
TO, see discussion in Sec. IV. 
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FIG. 11: (Color online) Elastic-elastic correlation length of 
the cross section at different values of A as a function of k 
for M = 1 and M = 10; the GOE case (connected circles), 
A = (pluses connected by a line), A ^ oo (squares), and 
A = 1/30 (crosses). The Weisskopf relation (I62|) is shown by 
dashed curve. 



The fact that the correlation length is not equal to the 
average width was recognized long ago, see Refe. JsiJ 
and [a]. However, the statements based on the equality 
(1391) still appear in the literature, see, for example, Ref. 
[4l|. The relation between the average resonance width 
and the transmission coefScient is given by the Moldauer- 
Simonius (MS) formula 



M\n{l-T) = -2tt 



(T) 
D 



(63) 



It can be seen from this expression and Eq. ([621) that the 
equality I = (F) is true only for small T. 



IX. CONCLUSIONS 

In conclusion, we have studied the statistics of cross 
sections for a fermion system coupled to open decay chan- 
nels. For the first time we carefully followed various sig- 
natures of the crossover from isolated to overlapping reso- 
nances in dependence on the strength of inter-particle in- 
teraction modelled here by the two-body random ensem- 
ble. The study was performed for the simplest Gaussian 
ensemble of decay amplitudes. Even in this limiting case, 
when these amplitudes were considered as uncorrelated 
with the intrinsic dynamics, we found significant depen- 
dence of reaction observables on the strength of intrinsic 
interaction. We expect this dependence to be amplified 
with realistic interplay of decay amplitudes and internal 
wave functions. Such studies should be performed in the 
future. 

A detailed comparison has been carried out of our 
results with standard predictions of statistical reaction 



theory. The average cross section was compared with 
the Hauser-Feshbach formula for a large number of chan- 
nels. In the inelastic case this description works quite well 
in the overlapping resonance regime for any interaction 
strength, while in the elastic case strong deviations have 
been found if the intrinsic motion is not fully chaotic. 

The study of Ericson fiuctuation theory shows that the 
assumption that the fiuctuations of the resonance widths 
become negligible for a large number of channels is wrong 
in the overlapping regime. We found that the fiuctua- 
tions of resonance widths increase with the coupling to 
the continuum, and we gave evidence that the relative 
fiuctuation of the width (the ratio of the variance to the 
square of the average width) diverges at k = 1 for any 
number of channels. This should imply that for any num- 
ber of channels the differences from the standard theory 
should increase as k increases. 

In order to study the relationship between the variance 
of the cross section and its average value, it is necessary 
to take into account the dependence of the average cross 
section on the intrinsic interaction strength A. Even when 
this is done, the standard prediction about the variance 
of the cross section was found to be a good approxima- 
tion only for a very large number of channels. For M 
between 10 and 20, where the Ericson prediction could 
be expected to be valid, consistent deviations have been 
demonstrated. In particular, the distribution of cross sec- 
tions shows that the probability of a large value of the 
cross section, mainly for the elastic case (or in the pres- 
ence of direct reactions), can be well below conventional 
predictions. 

Finally, we have shown that, in agreement with pre- 
vious studies, the correlation length differs from the av- 
erage width for any number of channels. On the other 
hand, the Weisskopf relation ([5^ that connects the corre- 
lation length of the cross section to the transmission coef- 
ficient, works, for a large number of channels, at any value 
of the intrinsic interaction strength A. In many situations 
we have seen that increase of A in fact suppresses the 
fiuctuations in the continuum. This can be understood 
qualitatively as a manifestation of many-body chaos that 
makes all internal states uniformly mixed. 

Our results can be applied to any many-fermion sys- 
tem coupled to the continuum of open decay channels. 
The natural applications first of all should cover neutron 
resonances in nuclei, where rich statistical material was 
accumulated but the transitional region from isolated to 
overlapped resonances was not studied in detail. The in- 
teresting applications of a similar approach to molecular 
electronics and electron tunneling spectroscopy can be 

Other open meso- 



found in the recent literature i^, 
scopic systems, for example, quantum dots and quan- 
tum wires, should be analyzed as well in the crossover 
region. One can expect very promising designated stud- 
ies for checking the statistical properties of resonances 
in such controllable experiments as those in microwave 
cavities and in acoustical chaos. Open boson systems in 
atomic traps also can be an interesting object of future 
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theoretical and experimental studies. 
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